Factors controlling the groundwater transport of U, Th, Ra, and Rn

A model for the groundwater transport of naturally occurring U, Th, Ra, and Rn nuclides in the U and Th decay series is discussed. The model developed here takes into account transport by advection and the physico-chemical processes of weathering, decay, -recoil, and sorption at the water-rock interface. It describes the evolution along a flowline of the activities of the U and Th decay series nuclides in groundwater. Simple sets of relationships governing the activities of the various species in solution are derived, and these can be used both to calculate effective retardation factors and to interpret groundwater data. For the activities of each nuclide, a general solution to the transport equation has been obtained, which shows that the activities reach a constant value after a distance xi, characteristic of each nuclide. Where xi is much longer than the aquifer length, (for U , U, and Th), the activities grow linearly with distance. Where xi is short compared to the aquifer length, (for Th,Th, Th, Ra, and Ra), the activities rapidly reach a constant or quasi-constant activity value. For Ra and Rn, the limiting activity is reached after 1 km. High U values (proportional to the ratio "234Th=w238U) can be obtained through high recoil fraction and/or low weathering rates. The activity ratios Th/Th, Ra/Ra and Ra/Ra have been considered in the cases where either weathering or recoil is the predominant process of input from the mineral grain. Typical values for weathering rates and recoil fractions for a sandy aquifer indicate that recoil is the dominant process for Th isotopic ratios in the water. Measured data for Ra isotope activity ratios indicate that recoil is the process generally controlling the Ra isotopic composition in water. Higher isotopic ratios can be explained by different desorption kinetics of Ra. However, the model does not provide an explanation for Ra/Ra and Ra/Ra activity ratios less than unity. From the model, the highest Rn emanation equals 2". This is in agreement with the hypothesis that Rn activity can be used as a first approximation for input by recoil (Krishnaswami et al 1982). However, high Rn emanation cannot be explained by production from the surface layer as formulated in the model. Other possibilities involve models including surface precipitation, where the surface layer is not in steady-state.